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Abstract. For known gravitational lens systems the red- 
shift distribution of the lenses is compared with theoreti- 
cal expectations for 10'* Friedmann-Lemaitre cosmological 
models, which more than cover the range of possible cases. 
The comparison is used for assigning a relative probabil- 
ity to each of the models. The entire procedure is repeated 
for different values of the inhomogeneity parameter rj and 
the limiting spectroscopic magnitude mam, which is im- 
portant for selection effects. The dependence on these 
two parameters is examined in more detail for the spe- 
cial cases Ao = and k = 0. 

Previous results that this method is a better probe 
for Ao than SIq are confirmed, but it appears that the 
low probability of models with large Aq values found using 
similar methods is due to a selection effect. 

The power of this method to discriminate between cos- 
mological models can of course be improved if more grav- 
itational lens systems are found. However, our numerical 
simulations indicate that a reasonable number of observed 
systems cannot deliver interesting constraints on the cos- 
mological parameters. 

Key words: gravitational lensing - cosmology: theory - 
cosmology: observations 



1. Introduction 

It has recently been suggested by many authors (see, 
for example, Fukugita et al. (1992) and references therein) 
that gravitational lensing statistics can provide a means 
of distinguishing between different cosmological models, 
most effectively concerning the value of the cosmologi- 
cal constant. This is fortunate, since most of the classi- 
cal methods for determining cosmological parameters are 
more sensitive to other quantities such as the density (Slo) 
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or deceleration (qo) parameters. It has even been sug- 
gested (Carroll et al. 1992) that gravitational lens statis- 
tics based on current observations already give the best 
upper limits on Aq for world models with k > 0, and are 
the most promising method of doing so for k = 0. 

Kochanek (1992) has suggested a method based not on 
the total number of lens systems but rather on the red- 
shift distribution of known lens systems characterised by 
observables such as redshift and image separation. Look- 
ing at a few different models, he concludes that fiat, A- 
dominated models are five to ten times less probable than 
more 'standard' models. The advantage of this method 
is that it is not plagued by normalisation difficulties as 
are most schemes involving the total number of lenses. 
The aim of this paper is to extend this method to arbi- 
trary Friedmann-Lemaitre cosmological models as well as 
to look at the infiuence of observational bias concerning 
the brightness of the lens. In addition, numerical simula- 
tions are used to estimate the usefulness of the method 
when more systems are available. 

It is important to note that the method described in 
this paper treats Aq and SIq as independent parameters, 
that is, they can in principle be determined simultane- 
ously. Also important is the fact that SIq is the global 
value, i.e., determined by the contribution of all com- 
ponents, regardless of degree of homogeneity and so on. 
This is because the cosmological parameters make them- 
selves felt through the cosmological model; most methods 
of determining SIq will miss any matter homogeneously 
distributed on a scale larger than that surveyed. 

2. Theory 

We make the 'standard assumptions' that the Universe 
can be described by the Robertson- Walker metric and 
that lens galaxies can be modelled as non-evolving singu- 
lar isothermal spheres (SIS). If one drops the first assump- 
tion, the cosmological parameters Aq, SIq and Hq lose their 
significance; the second assumption allows easy calcula- 
tion but, more importantly, is probably justified within 
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the attainable accuracy (see Krauss & White (1992) for 
a discussion). In order to have a well-defined statisti- 
cal quantity, which is based on the optical depth dr for 
'strong' lensing events, "'^ this discussion is limited to grav- 
itational lens systems with sources which are multiply im- 
aged (^ image separation) by isolated (^ negligible clus- 
ter infiuence) single galaxies and with known source and 
lens redshifts. An additional requirement is that the sys- 
tem must have been found without any biases concerning 
the redshift of the lens. (See Kochanek (1992) for a dis- 
cussion of these selection criteria.) 

Making use of the fact that the SIS produces a constant 
defiection angle, i.e., independent of the position of the 
source with respect to the optical axis (defined as passing 
through observer and lens), one can define the angular 
cross section ira"^ of a smgle lens for 'strong' lensing events 
(Turner et al. 1984): 

where v is the one-dimensional velocity dispersion of the 
lens galaxy, c the speed of light and Djs (-Ds) the angu- 
lar size distance between lens and source (observer and 
source). Following Kochanek (1992), one can arrive at an 
expression for the optical depth as follows. 

For a fixed mass and mass distribution (^ v), world 
model and Zg, the differential optical depth due to all 
lenses of a given mass as a function of z<j is of course 
proportional to the number of lenses of the given mass 
per z<j-interval and to the cross section for strong lens- 
ing events. In order to arrive at an expression for dr for 
a fixed image separation, one needs to know the relative 
number of lenses which, under the given circumstances, 
can produce this image separation. This can be done by 
using the Schechter luminosity function (Schechter 1976) 
as well as the Faber-Jackson and Tully-Fisher relations 
(Faber & Jackson 1976, Tully & Fisher 1977), which give 
the dependence of the velocity dispersion on the luminos- 
ity for elliptical and spiral galaxies, respectively. Bringing 
in the familiar parameters and dropping all terms which 
are concerned only with normalisation, one arrives, after 
some tedious but trivial calculations, at the expression 

dr ,2 « 7 / « A V^'^"^ 

= (1 -(- 2d) X 

dzd a* 2 \a* D^s ) 

where a* := 47r (") (t;* := v of an _L* galaxy), 7 is the 
Faber-Jackson/Tully-Fisher exponent, a the Schechter ex- 
ponent, _Dd the angular size distance between the observer 
and the lens and 

0(zd) := (1 + Zd)2(fio^d + 1 - Ao) + Ao. (3) 

^ See Schneider et al. (1992) for a clarification of the concept 
of optical depth in lensing. 



Equation (2) is independent of the Hubble constant since 
the dependences on iJo in the angular size distances 
and in the Faber-Jackson/Tully-Fisher relation cancel. 
In order to facilitate comparison with other authors, we 
have chosen the 'standard values' —1.1, 2.6, 4, 144 km/s 
and 276 km/s for the Schechter exponent, the Tully- 
Fisher exponent, the Faber-Jackson exponent, i'*spirai and 
t'*eiiipticai, respectively. (The value for t;*eiiipticai includes 
the factor (3/2)3 advocated by Turner et al. (1984) and 
so our elliptical galaxies correspond to the c = 2 models 
examined by Kochanek (1992).) 

The optical depth depends on the cosmological model 
through (5(zd) as well as through the angular size dis- 
tances, because of the fact that _Djj = Dij(zi, Zj] Aq, Slo, ??)• 
The infiuence of rj, which gives the fraction of homoge- 
neously distributed, as opposed to compact, matter is 
felt only in the calculation of the angular size distances, 
whereas the cosmological model in the narrower sense 
makes its infiuence felt here as well as through (5(zd). 

In general, there is no analytic expression for the Dij ; 
they can be obtained by the solution of a second-order 
differential equation. (See Kayser (1985) for the deriva- 
tion of the differential equation, also Linder (1988) for 
a more general formulation. For an equivalent derivation 
for Ao = see Schneider et al. (1992). Kayser et al. (1995) 
give a general discussion and an easy-to-use numerical im- 
plementation.) If one has an efficient method of calculating 
the angular size distances, it is easy to evaluate Eq. (2) for 
various world models described by the parameters Aq, SIq 
and rj. 

Worthy of note is the independence of Eq. (2) on the 
source luminosity function (which of course will generally 
itself depend on Zg as well), the relative numbers of galaxy 
types (the galaxy type for a particular lens is known) and 
the fraction of galaxies in clusters (the method looks only 
at field galaxies); these factors have to be taken into ac- 
count when doing statistics based on the total number of 
lenses. Also, Eq. (2) is insensitive to finer points of the 
mass model such as core radius and ellipticity (cf. Krauss 
k White (1992), Narayan k Wallington (1992)). The main 
idea is to compare the observed distribution of lens red- 
shifts with theoretical expectations for various world mod- 
els; the method is described in the next section. 

Equation 2 can take on appreciable values at interme- 
diate redshifts even though the lens galaxy would he too 
famt to he seen at the redshift m question (m > mam). 
In order to correct for this effect, we have calculated the 
redshift at which the lens galaxy would become too faint 
to have its redshift measured for the investigated cosmo- 
logical model and truncated dr at this point. (Details in 
the next section.) It is immediately obvious that failure to 
correct for the faintness of the lens galaxies will artificially 
exclude cosmological models with a high median redshift 
in Eq. (2), which might otherwise not be excluded. 
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Table 1. Gravitational lens systems used. For references see Refsdal & Surdej (1994) and references therein. Note that 9 
corresponds to the radius of the Einstein ring or half the image separation 
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3. Calculations 

The following gravitational lens systems meet our selec- 
tion criteria: 0142-100 (= UM 673 ), 0218+357, 1115+080 
(= Triple Quasar), 1131+0456, 1654+1346 and 3C324. 
(See Table 1 for observational data on these systems.) We 
considered the following ranges of values for the cosmo- 
logical parameters: 

-10 < Ao < +10 
< fio < 10 
< r] < 1 

which, of course, are much larger than contemporary 
knowledge demands. However, the history of cosmology 
shows that the knowledge of today is often out of fashion 
tomorrow, so that we prefer to develop an approach capa- 
ble of dealing with a wide range of cosmological models. 
Also important is the fact that it would be an additional, 
though by no means necessary, point in favour of the valid- 
ity of the method if it assigns the highest probability to a 
cosmological model within the presently accepted canoni- 
cal parameter space. 

We looked at 100 x 100 models in the Aq-SIq plane for 

r] = 0.0, 0.3, 0.5, 0.7, 1.0 for = 23.5 

and 

miim = 23.5, 24.5, cxd for rj = 0.5 



(Johnson R magnitudes). In addition, we looked at 100 x 
100 models in the ry-Slo and miim-Slo planes for the special 
cases of Ao = and k = 0. 

Before one can examine the relative probability of a 
given cosmological model, one must first see if it is com- 
patible with the observations. (Of course, this does not 
imply that the model is compatible with all observa- 
tions, merely with the ones necessary for this analysis: 
Zl, Zs, 0" and galaxy type.) One obvious restriction is that 
the largest source redshift Zg ^iax in the sample must be 
smaller than z^nax, the maximum redshift possible in the 
cosmological model in question. (See, e.g., Stabell & Refs- 
dal (1966) or Feige (1992) for a discussion of these cosmo- 
logical models.) Another restriction concerns the bright- 
ness of the lens. From the observables zi, Zg, 0" and galaxy 
type one can use Eq. (1) to calculate the velocity disper- 
sion V, transform this to an absolute luminosity using the 
Faber-Jackson or Tully-Fisher relation and then calculate 
the apparent magnitude for the given cosmological model 
(given by the angular size distance up to powers of (1 + zi) 
and A'-corrections^). If this calculated lens brightness for 



The apparent luminosity of the lens galaxy was calculated 
for the Johnson i?-band using the /('-corrections of Coleman, 
Wu & Weedman (1980) and applying a standard B — R cor- 
rection (since the 5-band Faber-Jackson and Tully-Fisher re- 
lations were used). Since these /('-corrections are based on dis- 
placement of standard spectra at 2 = which extend into the 
UV-band, they are given only up to 2 = 2.0, where evolution- 
ary effects would in any case have to be considered. However, 
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at least one lens is fainter than mam then the model is 
also incompatible with the observations. Since a realistic 
value for mam is at least a magnitude fainter than any 
miens value in Table 1, there is no danger that the actual 
cosmological model would be excluded by this restriction, 
even allowing for the uncertainty in calculating the lens 
brightness in the matter described. 

In these two cases we assigned the corresponding world 
model a probability of in our plots, indicated by white. 
Since the value doesn't occur otherwise, all white areas 
are due to these two restrictions. Otherwise, to measure 
the relative probability of a given cosmological model, we 
defined the quantity / as follows: 

< / := < 1' (4) 

Jo 

where zi is the observed lens redshift for a particular sys- 
tem. (z<j is used to denote the variable corresponding to 
lens redshift as opposed to the measured value for a partic- 
ular lens.) The distribution of the different / values (one 
for each lens system in the sample) in h equally-sized bins 
in the interval ]0,1[ gives the relative probability p of a 
given cosmological model, with 

1=1 

where rii is the number of systems in the i-th bin. This 
definition allows only a few discrete values, of course. The 
variable & is a free parameter; since the most information 
is obtained when h is equal to the number of systems, we 
adopt this value for h. 

Were the other observables the same for all systems, 
the redshift distribution should be given by Eq. (2); since 
this is not the case, the quantity / is defined, which al- 
lows one to compare the observed with the expected red- 
shifts for different observables and hence different relative 
probability curves (Eq. (2)) for each system. Our ansaiz 
is thus to expect that the /-values should be uniformly 
distributed for the correct cosmological model (barring in- 
trinsic scatter, of course). Simple combinatorics (the num- 
ber of ways to distribute a objects in h bins) and neglect- 
ing normalisation then leads to Eq. (5). For a given world 
model, the 6 /-values (one for each gravitational lens sys- 
tem used) are calculated, and these values are used to 
determine the relative probability via Eq. (5). 

4. Results and discussion 

Figure 1 shows the relative probability of several cosmo- 
logical models as given by Eq. (5). The grey scale in all 

in most cases the galaxy becomes too faint at modest redshifts, 
so the assumption of no evolution made throughout this paper 
is probably justified. 



Table 2. Values of the relative probability in the plots in Fig. 1 
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mum values of each individual plot. The resolution is 0.2 
in Ao and 0.1 in SIq, thus giving 10.000 different models. 
The scale is at the right. (Because the relative probabil- 
ity can only take on a few discrete values, contour levels 
aren't very useful as indicators of the relative probability, 
since they would merely indicate the boundaries between 
areas of constant probability. In order to indicate the val- 
ues directly. Fig. 1 plots the relative probability on a grey 
scale. Although in themselves not important, the inter- 
ested reader can read off the probabilities directly in the 
legend, where the discrete values which actually appear in 
the plots have been marked. In addition. Table 2 gives the 
values of the relative probability which occur in each plot.) 
All plots except (9), (10), (12) and (13) are in the Aq-Qo 
plane. Plot (8) gives some orientation in this plane. The 
diagonal line from upper left to lower centre corresponds 
to A; = 0; the six curves are, left to right, for htg = 4, 5, 
6, 8 and 10 x lO^a (h := Hq ■ lOO-^km"! • s • Mpc) as well 
as the border to the so called bounce models, i.e. models 
with no big bang and thus a maximum redshift z^nax', the 
line on the right corresponds to qo = —5. The values of 
the fixed parameters rj and mum are indicated on each 
plot. Although some of the area in this plane is definitely 
excluded^ by simple arguments, a probability has been 
computed for each world model, in keeping with the sec- 
ond of our motivations mentioned above. In the white area 
at the right the probability is because these world models 
have a maximum redshift lower than the highest redshift 
in Table 1; in plots (1), (2), (3) and (7) there is an addi- 
tional white area (p = 0) separated from the first one by a 
strip where p > due to the fact that at least one lens is 
fainter at its observed redshift than mum. (The brightness 
of the lens was not used as an additional constraint due 
mainly to the fact that the computed brightnesses are only 
correct to about a magnitude or so (Kochanek 1992).) 

Even though not every point in the plane, i.e. every combi- 
nation of Ao and Qo, corresponds to a world model which can- 
not be excluded by simple arguments, nevertheless the ranges 
of the individual parameters are allowed assuming the lowest 
realistic values for Ho and the age of the universe. 
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(9) k = 0, miinj = 23.5 
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Fig. 1. Results for the systems in Table 1. The relative probability for different cosmological models is plotted linearly on 
the scale shown at the right, where the discrete values which appear in the plots are marked. Fixed parameters are indicated 
above each plot. For comparison, the results for the systems used in Kochanek (f992) are included in plot (6). Plot (8) is for 
orientation 
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Looking at plots (1), (2), (3), (4), (5), (7) and (11) the 
general impression is that, apart from some small-scale 
structure at the lower right of the p > area which ap- 
pears to be significant, there is not much structure. Since 
there are only discrete values of p, there are some dis- 
continuities. Apart from this large-scale-structure, how- 
ever, there is not much change from world model to world 
model. The curves along which the probability is constant 
are more vertical than horizontal, suggesting that this 
method better probes Aq than SIq. Although not of any 
use statistically, for each plot in Fig. 1 the Kolmogorov- 
Smirnov probability (see Sect. 5) was also computed. In 
the interest of brevity, we have not included these plots 
here; however, the structure is qualitatively similar (with- 
out the sharp discontinuities, of course) and also hints 
at the significance of the small-scale structure mentioned 
above. 

Plots (1), (2), (3), (4) and (5) show the efl'ect of vary- 
ing rj. As one can see, the general structure doesn't change 
much, especially if one concentrates on the area near the 
roughly vertical thin strip (which would still be allowed 
assuming higher values for Hq and/or the age of the uni- 
verse) for relatively low SIq values (SIq < 4, say). This is 
compatible with the result of Fukugita et al. (1992) which 
indicates only a weak dependence of the lensing cross 
section on rj. rj enters into the calculation only in the 
computation of the angular size distances, whereas Aq 
and Slo also enter the calculation through the function 
Q(zd) as defined in Eq. (3). Also, the particular combina- 
tion of the angular size distances involved in computing 
the lensing cross section, (_D<j-Dds)/-Ds (see, e.g., Fukugi- 
ta et al. (1992) or Kochanek (1992); in Eq. (2) the explicit 
dependence on _D<j has been lost in the necessary variable 
transformation and integration (cf. Kochanek (1992))), is 
relatively insensitive to rj. This is not true, for exam- 
ple, for the combination _Dds/(-Dd-Ds), which is impor- 
tant for computing Hq from the time delay between the 
different images of a multiply imaged source (cf. Fukugi- 
ta et al. (1992)). 

Plots (3), (7) and (11) show the effect of varying mam. 
Of course, the additional white area described above be- 
comes smaller as mam becomes fainter, disappearing for 
miim oo (this condition is sufficient but not necessary 
for every lens galaxy to be brighter than mum for the world 
models examined). The main differences, however, occur 
in the small-scale structure at the lower right of the p > 
area: the fainter mum becomes, the less probable the mod- 
els near the border of this area appear. For example, if one 
compares the models near the Einstein-de Sitter model 
(Ao = and SIq = 1) with the models near the de Sitter 
model (Ao = 1 and SIq = 0), two models which have been 
examined rather extensively in the literature on lensing 
statistics, especially in the direction 'perpendicular' to the 
curves dividing different probabilities, then one sees that 
for miini — 23.5 and mi^m — 24.5 those near the de Sitter 
models have a roughly equal but slightly higher relative 



probability than those near the Einstein-de Sitter model; 
only for mum = cxd is the situation reversed, those near 
the Einstein-de Sitter model having a clearly higher prob- 
ability. 

This is easy to understand, since it is these models 
near the de Sitter model which have a maximum in dr 
at relatively large redshifts (cf. Kochanek (1992)); for re- 
alistic values of miim, one cannot see the lens galaxies at 
these redshifts. If one uses a realistic value for mum, one 
compares the redshift distributions for the different world 
models (given by Eq. (2)) at small redshifts, where they 
don't differ very much. If one takes mum = cxd, implying 
that one could measure the redshifts of the lenses at all 
redshifts, whatever their brightness, then it appears that 
models with a large median redshift in Eq. (2), such as 
those on the right hand border of the p > area includ- 
ing the models with a large cosmological constant exam- 
ined by Kochanek (1992), are improbable. However, this is 
merely a selection effect. The relatively low lens redshifts 
m Table 1 don't mean that world models with a large me- 
dian lens redshift are improbable; it means that we can't see 
the lenses at these redshifts. Comparing probability distri- 
butions assuming that we could artificially excludes these 
models m preference to models like the Emstem-de Sitter 
model with a small median lens redshift. Introducing mum 
makes the situation more realistic, but means comparing 
the distributions at small redshifts. Thus, the power to 
discriminate between various world models is reduced. 

For realistic values of mum, such as those in plots (3) 
and (7), the exact value of mum isn't very important; what 
does make a difference is assuming a value which is much 
too faint. This is easy to understand, since, near the red- 
shift where a lens galaxy typically becomes fainter than a 
realistic value for mum, the function m(z<j) is rather steep. 
This means that a change in mum by a magnitude or two 
corresponds to a relatively small change in z<j and thus to 
a correspondingly small change in the area under the dr 
curve up to this z<j value; thus, there is little infiuence on 
the value of / as defined in Eq. (4). This is illustrated in 
Fig. 2. 

For comparison, we have also tested the method on 
the systems used by Kochanek (1992), using mum = cxd 
und = 1, both of which he implicitly assumes.'* The re- 
sults are in plot (6) where the relative probabilities are ^ , 
i and 1 and comparing the various models examined by 
Kochanek confirm his conclusions. For instance, the rela- 
tive probabilities of the Einstein-de Sitter and de Sitter 
model are 1 and confirming his result that fiat, A- 
dominated models are 5-10 times less probable than stan- 
dard ones. (However, taking mum into account and/or us- 



Of course, when one considers finite values for mum, one 
cannot include systems with lens redshifts which have been 
determined by means other than measured emission redshifts, 
such as absorption lines (which assumes that the lens is also 
the absorber). 




Fig. 2. The relative differential opticai depth (thin curve) and 
the cafcufated fens brightness m (thicfc curve) as functions of 
Zd- The worfd modef is the de Sitter modef (Ao = f .0 Qo = 0.0; 
the vafue of rj doesn't matter since there is no matter) and 
the observabies are those for the gravitational lens system 
0142 — 100 (see Table 1). The ordinate gives the magnitude 
in Johnson R 

ing only directly measured lens redsliifts would produce 
quite different results, as discussed above.) This plot ar- 
tificially indicates a low probability for models near the 
de Sitter model for the same reasons as those discussed in 
connection with plot (11). 

Plots (9) and (10) examine the influence of ry and miim, 
respectively, for the special case k = 0. Plots (12) and 
(13) examine the influence of rj and mam, respectively, for 
the special case Aq = 0. In plots (9) and (12) one can 
easily see the weak dependence on rj, especially for small 
values of SIq. In plots (10) and (13) one can see the even 
weaker dependence on mam in this range. There are hints 
toward fainter magnitudes of a declining probability for 
small Slo values, as discussed above. The white area at 
the left in plots (10) and (13) is due to the fact that at 
least one lens is fainter at its observed redshift than the 
corresponding mum value. (Since these mum values are 
unrealistically small, these models are not incompatible 
with the observations.) 

At first glance, the 'oscillations' in the relative prob- 
ability might appear somewhat puzzling. According to 
Eq. (5), a higher probability is obtained for a more regular 
distribution. Since Eq. (5) only allows discrete values for 
the relative probability, a 'jump' occurs when the num- 
ber of systems in a certain bin changes (unless offset by 
a corresponding change in another bin). This can be seen 



Fig. 3. 'Oscillations' in the relative probability. The thin hori- 
zontal lines are the bin boundaries. The thick curves show, from 
top to bottom, the /-values given by Eq. (4) for the gravita- 
tional lens systems 1131-|-0456, 3C 324, 0218-1-357, 0142-100, 
1115 -|- 080 and 1654 -|- 1346, respectively. Since two systems 
have such large lens redshifts that dr is practically = at the 
corresponding redshift, their /-values (in these world models) 
are practically = 1 and cannot be distinguished from each other 
or from the bin boundary at / = 1 in the resolution of the plot. 
Table 1 shows that these are the two highest lens redshifts in 
the sample. Beneath the lines and curves, shown in the same 
way as in Fig. 1, is the corresponding relative probability 

in Fig. 3, where for demonstration purposes the /-values 
given by Eq. (4) for each gravitational lens system used 
are plotted as functions of Aq for k = and rj = 0.3. 
(That is, for the world models along the k = line in 
plot (4) in Fig (1).) It can be seen that, although — due to 
the definition — the relative probability changes by notice- 
able amounts between a few discrete values, nevertheless 
the /-values themselves are smoothly varying functions of 
the world model. Figure 3 also makes the following general 
conclusions clear in this specific example. 

— The fact that a couple of systems have /-values which 
are practically = 1 limits the maximum probability, 
since these are always in the same bin. 

— 'Oscillations' between different probabilities have no 
physical significance, but rather are merely artifacts 
of the particular lens redshifts. On the other hand, 
extremely low probabilities, e.g. all six systems in the 
same bin (p = ^^g), would be more indicative of a 
low-probability cosmological model. 

— Apart from the oscillations, a trend (the /-values in- 
creasing to the left in Fig. 1) can be seen which would 
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indicate a probability low enough to reject the cor- 
responding cosmological models, were the parameter 
space examined larger. Thus, the method might be 
able to exclude 'extreme' cosmological models. 

The interested reader can use Fig. 3 together with Eq. (5) 
to see how the relative probability is arrived at. 

The statistical significance of all results in this section 
is discussed in the next section. 

5. Numerical simulations 

Numerical simulations were done for mam = 23.5 and 
for rj = I. This value of mam = 23.5 is the most realis- 
tic based on the present state of observations and using 
only one value for mum = 23.5 as well as for rj is justi- 
fied based on the weak dependence of the results on these 
parameters, as discussed in Sect. 4. 

The observables 6" (the radius of the Einstein ring or 
half the image separation corresponding to the diameter 
of the Einstein ring), Zg and galaxy type were chosen ran- 
domly from an interval roughly corresponding to the ob- 
served range of values in order to produce synthetic data 
comparable to real observations. It is important to note 
that neither the exact range nor the shape of the distribu- 
tion matters, since the method looks at the redshift dis- 
tribution of the lenses with the other parameters fixed by 
observation. For convenience, a fiat distribution was cho- 
sen for each of the observables. For a given cosmological 
model, the corresponding lens redshift zi for each system 
was calculated from the observables and a randomly gen- 
erated / through (numerical) inversion of Eq. (4). This 
catalog was then used to determine a relative probabil- 
ity for each of the points in the Aq-SIq plane in the same 
manner as for the real systems. 

With the probability given by Eq. (5), based on sim- 
ple combinatorics, one cannot know to what degree the 
values for each cosmological model are infiuenced by sta- 
tistical fiuctuations. However, with such a small number 
of systems, there is really no other method of computing a 
relative probability. We expect that the distribution of the 
/-values, barring statistical fiuctuations, should be fiat for 
the correct cosmological model. So we need a test to com- 
pare this distribution with a fiat probability distribution. 
The Kolmogorov-Smirnov (K-S) test is of course a well un- 
derstood method for testing if two distributions are statis- 
tically significantly different. (See, e.g.. Press et.al. (1992) 
for a general discussion and definition of the K-S probabil- 
ity.) However, this test can only be used for distributions 
with more than 20 data points. For purposes of com- 
parison, for the systems in Table 1 not only was the prob- 
ability defined in Eq. (5) computed (shown in Fig. 1) but 
also the K-S probability. The K-S probability should of 
course not be taken seriously if there are too few systems. 

We have done simulations for a variety of world mod- 
els and also for numbers of systems between 20 and 50. In 
the interest of brevity, we present only one plot. Figure 4 
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Fig. 4. Results based on a catalogue of 50 simulated systems. 
(Note the different scale on the axes.) The cosmological model 
used to generate the lens redshifts is the homogeneous Ein- 
stein-de Sitter model {rj = 1.0, Ao = 0.0, Qo = 1.0). The grey 
scale is as in Fig. 1 

shows the results derived from a catalogue of simulated 
gravitational lens systems. Since, even with 50 systems, 
no area can be excluded based on the K-S probability - 
the white area has p = due to the fact that at least 
one lens would be fainter than mum in these world mod- 
els, as discussed above - we conclude that, although one 
can qualitatively understand the physics which at least in 
part is responsible for the results presented in Fig. 1, the 
actual relative probabilities are more indicative of intrin- 
sic scatter in the redshifts of the lenses than a hint of the 
correct cosmological model. 

6. Summary and conclusions 

In this paper we have extended the method originally pro- 
posed by Kochanek (1992) for using the redshift distri- 
bution of gravitational lenses to learn something about 
the cosmological model. This method is particularly at- 
tractive since it avoids the normalisation difficulties nor- 
mally associated with lensing statistics. First, we made 
use of the equation derived by Kayser (1985; see also 
Kayser et al. 1995) to be able to examine cosmological 
models described by arbitrary values of Aq, SIq and rj. Sec- 
ond, we looked at the infiuence of observational bias by 
introducing the parameter mum. (This means that we can- 
not include systems where the lens redshift has been esti- 
mated by some other means than an emission spectrum.) 
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Third, we used a more quantitative statistic to look at 
relative probabilities. 

Since the potential reward from using the Kochanek 
formalism (setting limits on cosmological parameters bet- 
ter than with other methods using only observable quanti- 
ties and standard assumptions whose validity in this con- 
text is undisputed) seemed large, our basic idea was to 
extend this formalism to enable it to look at a larger num- 
ber of cosmological models (arbitrary values for SIq and 
Ao as well as rj) while correcting an obvious limitation 
(selection effects due to the brightness of the lens) and 
using a variation of his statistic (Kochanek (1992) basi- 
cally uses our statistic with two bins) which allows more 
information about the distribution of the redshift values 
to be used. That is, we intended to follow the formalism 
of Kochanek (1992) as closely as possible. Unfortunately, 
the fact that the relative optical depth given by Eq. (2) is 
appreciably different only for those cosmological models 
in which this difference cannot be seen due to the selec- 
tion effect renders the technique less useful than we had 
hoped. 

We saw that little extra uncertainty in the derived val- 
ues for Ao and SIq is introduced by letting ry be a free pa- 
rameter. The same is true of mam with the exception that 
values which are unrealistically faint distort the results. 
A comparison with the results of Kochanek (1992), con- 
firmed with our formalism, also show the consequences of 
neglecting mam. 

An interesting result is the degeneracy of the derived 
relative probability along curves roughly parallel to curves 
of constant world age in the Aq-SIq plane. This indicates 
that the method is more sensitive to Aq than to SIq, but 
also shows that demonstrating the consistency of a given 
cosmological model with the observations using this sta- 
tistical method also implies consistency with a large range 
of other cosmological models. 

The dramatic difference caused by not neglecting mum 
casts doubt on the degree to which present observations, 
based only on the redshift distribution,® are able to rule 
out certain cosmological models; in particular, fiat models 
with a large cosmological constant, having a high median 



More information is available in theory by looking at not 
only the redshift distribution, i.e., the shape of the curve, but 
also the number of lenses, i.e., the area under the curve. This, 
however, introduces additional uncertainties due to normalisa- 
tion. Additionally, one could consider the completeness of the 
sample, as pointed out to the authors by the referees. For ex- 
ample, the number of lens systems with unmeasured redshifts 
could in principle be used to exclude cosmological models in 
which the redshifts could have been measured, or the other 
way around. This requires information regarding the reasons 
as to why the redshifts haven't been measured — the effects of 
the cosmology should be separated from the current observa- 
tional stand, which requires a detailed analysis of the observa- 
tional literature, or a separate observational programme, both 
of which are beyond the scope of this paper. 



expected lens redshift, become more probable through in- 
troducing miim. 

Nevertheless, plots (6) and (11) gives an idea of what 
could be done, if one were able to measure the redshifts of 
the faintest lens galaxies. For a given image separation, the 
calculated brightness of the lens galaxy has a minimum at 
some intermediate redshift; this is typically at about 30™ 
in R, so that larger telescopes and advances in image pro- 
cessing will probably be able to make some progress on 
this front in the next several years. If one were able to 
measure the lens redshift at the minimum brightness, this 
would have the side-effect of eradicating the dependence 
on miim. On the other hand, probably more would be lost 
than gained, because it would no longer be possible to ne- 
glect evolutionary effects. For this reason, one could sup- 
pose that more progress in the immediate future (barring 
a revolution in the understanding of evolutionary effects) 
will probably come from increasing the number of usable 
systems (through the discovery of more systems and/or 
through measuring more redshifts in known systems) than 
from pushing mum to fainter values. 

The results of using the method on synthetic data, 
however, cast doubt on the statistical significance of our 
results and on the hope of using this method to exclude 
certain cosmological models, especially those which cannot 
be excluded by other tests. Thus, being conservative in our 
appraisal of what the statistics of redshifts of gravitational 
lenses can tell us, we conclude that at present and in the 
foreseeable future this method will probably not give us 
any useful information. 
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